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ABSTRACT 



Aims. The aim of this work is to search for radio signals in the quiescent phase of accreting millisecond X-ray pulsars, in this way 
giving an ultimate proof of the recycling model, thereby unambiguously establishing that accreting millisecond X-ray pulsars are the 
progenitors of radio millisecond pulsars. 

Methods. To overcome the possible free-free absorption caused by matter surrounding accreting millisecond X-ray pulsars in their 
quiescence phase, we performed the observations at high frequencies. Making use of particularly precise orbital and spin parameters 
obtained from X-ray observations, we carried out a deep search for radio-pulsed emission from the accreting millisecond X-ray pulsar 
XTE J0929-3 14 in three steps, correcting for the effect of the dispersion due to the interstellar medium, eliminating the orbital motions 
effects, and finally folding the time series. 

Results. No radio pulsation is present in the analyzed data down to a limit of 68 jxJy at 6.4 GHz and 26 yujy at 8.5 GHz. 
Conclusions. We discuss several mechanisms that could prevent the detection, concluding that beaming factor and intrinsic low 
luminosity are the most likely explanations. 

Key words. Pulsar: general - Pulsar: individual (XTE J0929-314), neutron star. X-ray binary 



1. Introduction 

The neutron star low-mass X-ray binaries (NS-LMXB) are sys- 
tems containing a neutron star (NS) believed to have a weak 
magnetic field (B < 10"^ G) and accreting matter from a low- 
mass (M < 1 M 0) companion. A special subgroup of NS- 
LMXBs is that of the neutron star soft X-ray transients (NS- 
SXT; see e.g. White et al. 1995). These systems are usually 
found in a quiescent state with luminosities in the range 10^' - 
10^^ erg ■ s-'. On occasion they exhibit outbursts with peak lu- 



minosities in the 0.5-10 keV range between 10 and 10 erg ■ 
s"', during which their behaviour closely resembles that of the 
' persistent NS-LMXBs (Campana et al. 1998). According to the 
recycling scenario (Alpar et al. 1982; Bhattacharya & van den 
Heuvel 1991), the NS-LMXBs are the progenitors of the mil- 
lisecond radio pulsars (MSPs), reaccelerated by mass and angu- 
lar momentum transfer from the companion star. 

Although widely accepted, there is no direct evidence to con- 
firm this scenario yet. However, in 1998, the idea that NSs in 
NS-LMXBs are spinning at millisecond periods was spectacu- 
larly demonstrated by the discovery of coherent X-ray pulsation 
at 2.5 ms in SAX J1808.4-3658, an NS-SXT with an orbital pe- 
riod Porb ~ 2 hr (Wijnands & van der KUs 1998; Chakrabarty & 
Morgan 1998). For almost four years, SAX J1808.4-3658 has 
been considered a rare object in which some peculiarity of the 
system (e.g. its inclination) allowed detection of the NS spin; 
however, the situation is now dramatically changed as nine other 
NS-SXTs in outburst have been discovered in which coherent 
X-ray pulsations in the millisecond range have been found (e.g. 
Wijnands et al. 2006; Ki'imm et al. 2007; Casella et al. 2008; 



Altamirano et al. 2008). We are, therefore, facing a new class of 
astronomical objects (dubbed accreting millisecond X-ray pul- 
sars: AMXPs) that may constitute the bridge between accretion- 
powered rapidly-rotating and rotation-powered NSs; in particu- 
lar, the detection of radio pulsations from these sources during 
quiescence would be the ultimate proof of the validity of the re- 
cycling model. 

During the mass transfer, the plasma from the companion 
settles into an accretion disk, whose inner rim is located at a 
distance rn, from the NS centre, with r^ the so-called magneto- 
spheric radius, which is equal to a fraction < 1 of the Alfven 
radius, at which the ram pressure of the infalling matter from the 
accretion disk balances the pressure of the NS magnetic field : 



1.0 X 10" 



cm. 



(1) 



Here yU26 is the magnetic moment in units of 10^^ G ■ cm-*. Re 
is the NS radius in units of 10^ cm, rhe the Eddington accretion 
rateinMo- yr"' (forRg = l,mE = 1.5xlO"^M0- yr"' and scales 
with the radius of the compact object), and mi the NS mass in 
solar masses. 

For the radio emission mechanism to switch on in a rotating 
magnetized NS, the space surrounding the NS must be free of 
matter up to the light cylinder radius (at which the speed of the 
material co-rotating with the NS would be equal to the speed of 
light): 



ric =5x10*' P- 



cm 



(2) 



where P_3 is the pulse period in milliseconds. 
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During an outburst, the occurrence, in some cases, of type 
I bursts and the observation of coherent X-ray pulsation indi- 
cate that at the origin of the observed luminosity is the accretion 
mechanism onto the NS surface. As a consequence coherent ra- 
dio emission cannot occur in this phase. In the quiescent state, 
the 0.5-10 keV luminosity is detected at levels ranging between 
lO-'''^ and 10-''* erg ■ s"' . Accretion of matter at a lower rate was 
originally proposed to account for this lower luminosity emis- 
sion. However detailed studies of the X-ray spectrum in quies- 
cence (Rutledge et al. 2001) and of the thermal relaxation of 
the NS crust during this phase (Colpi et al. 2001) suggest that 
the cooling of the periodically warmed up NS surface (Brown, 
Bildsten, & Rutledge 1998) is a viable explanation for the bulk 
of the luminosity in quiescence. If this is the correct interpreta- 
tion, there should be no mass accreted onto the neutron star sur- 
face during the X-ray quiescent phase, hence rfi ~ 0, and thus, 
plausibly. 



> ric 



(3) 



The timescale for the expansion of the magnetospheric radius 
beyond the light cylinder radius, in response to a sudden drop in 
the mass transfer rate, is much shorter (Burderi et al. 2001) than 
the typical duration (approximately years) of a phase of quies- 
cence in an NS-SXT, allowing the radio pulsar to switch on in 
principle. 

1.1. Absorption effect of the matter surrounding tine system 

In 2001 Burderi et al. proposed a model able to explain the fail- 
ure of previous and subsequent searches (Burgay et al. 2003 and 
references therein) for pulsed radio emission from quiescent NS- 
SXTs. As mentioned above, a temporary significant reduction of 
the mass-transfer rate may cause the magnetospheric radius to 
overcome ric, thus allowing a radio pulsar to switch on. In some 
cases, even if the secular mass-transfer rate is restored, the ac- 
cretion of matter onto the NS can be inhibited because the radi- 
ation pressure from the rotating dipole may be capable of eject- 
ing most of the matter overflowing from the companion out of 
the system. This phase has been defined as radio-ejection. One 
of the strongest predictions of this model is the presence, during 
the radio-ejection phase (hence during X-ray quiescence), of a 
strong wind of matter emanating from the companion star swept 
away by the radiation pressure of the pulsar (see Di Salvo et al. 
2008a for a discussion of a possible secular evolution of these 
systems). This matter, as well as that residual from a previous 
outburst, can play a role in absorbing the radio signal. Following 
Burderi et al. (2001) and Burgay et al. (2003), we can estimate 
the optical depth, Tff, at various radio frequencies due to matter 
engulfing an NS-SXT in its quiescent phase: 



Tff 



,n3 r'rfi_io(X + 0.5Y)2F(mi,m2)Ga 
1.6x10 X 



mf Tf Pf v2(l+m,/m2)^/3 



(4) 



where rh_io is the mass transfer rate in outburst in units of 
lO"'" Mq- y"', X and Y are the mass fraction of hydrogen and 
helium respectively, y is the fraction of ionized hydrogen, T4 the 
temperature of the outflowing matter in units of 10^ K, vg the 
frequency of the radio emission in units of 10^ Hz, Pj, is the or- 
bital period in hours, mi and m2 are the masses of the NS and its 
companion in solar masses, F (mi, mi) = 1 - 0.462 m2 / (mi + 
mj), Ga = Ga (T4, vg, Z) = 1.00 H- 0.48 (logT4 - log Z) - 0.25 
log vg takes into account the dependencies of the Gaunt factor 
from the temperature, the atomic number Z, and the frequency. 



Table 1. Orbital and spin parameters for XTE J0929-314 (Di 
Salvo et al. in prep.; see also Di Salvo et al. 2008b). The errors on 
the last quoted digit(s) are intended to be at 1 cr and are reported 
in parentheses. 



Parameter 


Value 


RAJ 


09" 29"" 20M9 


DECJ 


-31° 23' 3".2 


Orbital period, Po^b (s) 


2614.748(3) 


Projected semi-major axis, a- sin i (It-ms) 


5.988(10) 


Eccentricity, e 


< 0.007 


Spin period", Ps (s) 


0.0054023317856(4) 


Mean spin period derivative, Ps (s • s"') 


1.63(12) X lO^'* 


Ascending node passage, To (MJD) 


52405.48676(1) 



' The value reported is referred at epoch 52396.5 MJD. 



One example of a system in the radio-ejection phase is PSR 
J1740-5340, an eclipsing MSP, with a spin period of 3.65 ms 
located in the globular cluster NGC 6397 (D' Amico et al. 2001). 
For this source Tg = 2.8 at vg = 1.4 (with X = 0.7, Y = 0.3, 
y = 1, T4 = 1, Ga = 0.94 for Z = 1.1). This shows that the 
effect of free-free absorption at 1.4 GHz is important, but not se- 
vere, for this system. The radio signal from this PSR is indeed 
visible at this frequency although it is eclipsed and distorted at 
many orbital phases (D'Amico et al. 2001). X-ray millisecond 
pulsars are very similar to PSR J1740-5340, the main differ- 
ence being the orbital period, which is ~ 32 hr in the case of 
PSR J1740-5340 and typically ranges from 40 min to ~ 4 hr 
in the case of the X-ray millisecond pulsars. The difference in 
the orbital period results in a very different optical depth for the 
free-free absorption caused by the matter present around these 
systems. But the dependence of Tff on the square inverse of the 
radio frequency implies that, at higher frequencies, the effect of 
the surrounding matter will be less important. For instance, ap- 
plying Eq. Hto one of the tighter orbits AMXPs, XTE J0929- 
314, at 1.4 GHz, at 6.4 GHz and at 8.5 GHz, with parameters 
m_io ~ 2.9, mi ~ 1.4 Mq, mj ~ 0.02 Mq, X = 0.7, Y = 0.3, y ~ 
1, T4 ~ 1, Ph = 0.73 (Galloway et al. 2002) leads to Tff (1.4 GHz) 
w 5, Tff (6.5 GHz) X 0.2, Tff (8.5 GHz) ^Q.l. Adopting these two 
observing frequencies, the problem of the absorption of the radi- 
ation is totally overcome. Prompted by these considerations, we 
have undertaken a programme to search for millisecond pulsa- 
tions in NS-SXT XTE J0929-314 at 6.5 and 8.5 GHz. 

The observations and the method of data analysis are pre- 
sented in Sec.|2] and the results are reported in ^and discussed 
in an 



2. Observations and data analysis 

Radio observations of the millisecond pulsar XTE J0929-314 
were made with the Parkes 64 m radio telescope. Three data se- 
ries were taken on 2003 December 19-21 with a bandwidth of 
576 MHz. The first was at a central radio frequency of 6410.5 
MHz, the other two at a central radio frequency of 8453.5 MHz. 

The collected signal for each polarization was split into 192 
3 MHz channels using an analogue filterbank, in order to mini- 
mize the pulse broadening caused by dispersion in the interstel- 
lar medium (ISM). The outputs from each channel were one bit 
digitized every 100 i-is. The resulting time series were stored on 
Digital Linear Tapes (DLT) for off-line analysis. The observa- 
tions lasted 7.5 hr each, corresponding to 2^^ samples. We then 
rebinned the data to obtain 2^^ samples in order to reduce the 
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computational time; therefore, the effective time resolution of 
the analysed data is 200 ;us. 

The data analysis methodology was chosen on the basis of 
the precise knowledge of the orbital and spin parameters for 
XTE J0929-314, obtained from the X observation. The original 
ephemerides were presented by Galloway et al. (2002), while we 
used those reported in Table|2] which have been recently refined 
by Di Salvo et al. (in prep.; see also Di Salvo et al. 2008b). 

Off-line analysis was made by means of a software suite 
that, in the first stage, aimed to reduce the dispersion effects on 
the signal due to the interstellar medium. The data series were 
dedispersed according to 72 trial dispersion-measure (DM) val- 
ues ranging from 5.51 to 396.74 pc- cm"^ for the data series 
at 6410.5 MHz and 33 trial DM values ranging from 12.64 to 
416.99 pc- cm"^ for the data series at 8453.5 MHz. The max- 
imum DM was chosen to be ~ 4 times the nominal DM value 
obtained using either the Taylor & Cordes (1993) or the Cordes 
& Lazio (2001) models for the distribution of free electrons in 
the ISM, in order to account for the errors in the estimated dis- 
tance of the source, for the uncertainties in the ISM model, and 
for the presence of local matter surrounding the system, hence 
increasing the local DM (estimated to be at most ~ 100 pc/cm-' 
for this source; Burgay et al. 2003). 

The successive step is to deorbit and barycentre the time se- 
ries, i.e. to eliminate the effects of the orbital motion of the NS 
in the binary system and the earth in the solar system. This was 
obtained by resampling the time series in order to mimick time 
series collected from a telescope located in the barycentre of the 
solar system and looking at the source as if it were located at the 
barycentre of the pulsar orbit. 

In doing this exercise it is important to estimate the effect on 
the putative pulsar signal due to the uncertainties in the adopted 
ephemeris for the system. Therefore we simulated dedispersed 
time series containing a periodic signal with the spin and orbital 
characteristics of XTE J0929-314, as derived by Di Salvo et al. 
(Table |2]i. We then deorbited and barycentred these time series 
using the upper and the lower limits of the parameters in Table 
|2] adopting Icr errors. 

It turned out that even a variation in all the parameters (but 
the orbital period, see later) at Icr level does not affect the de- 
tectability of the pulsations, since that produced a maximum 
broadening of the pulse much smaller than 0. 1 in pulsar phase. 
The only exception is for the orbital period: Icr variation in Poit, 
propagated over the ~ 19500 orbits occurred between X-ray and 
radio observations, produces a broadening of 0.4 in pulsar phase 
for our observation. Therefore, to obtain a maximum broadening 
of the pulse of 0. 1 in phase, we corrected the time series using 
8 trial values of Pmb (4 above and 4 below the nominal value), 
covering the 1 cr uncertainty range. 

The third step in the procedure was to fold the (already de- 
orbited and barycentred) time series according to the rotational 
parameters reported in Table |2l The trial values of the spin pe- 
riod used for the folding, were chosen taking into account that 
the nominal value of Ps dates from the X ray observations of 
May 2002, while our observations occurred 19 months later. 
Therefore we explored a range of spin period APtot = Pmax- Pmin, 
with P„,a^ = Ps -H 5Ps + AT (Ps + dPs) and Pn,i„ = Ps - (5Ps, where 
Ps and Ps are the values reported in Table |2] 6Ps and 6Ps their 1 
cr errors, and AT is the time between the observations in X-ray 
and radio bands. As a first guess we used the outburst value of 
Ps for this calculation, resulting a safe upper limit for the dipolar 
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spin-dowrOl as shown below. For a safer investigation, we also 
checked on spin period values lower than our previous estimate 
Pmin and for the double of the nominal value, 2xPs. 

Since the spin period derivative measured in outburst is usu- 
ally different from that in the quiescent phase, we checked the 
plausibility of the adopted trial period interval through an esti- 
mate of the surface magnetic field Bs. This estimate can be ob- 
tained evaluating the magnetic torque acting on the NS through 
the formula (Rappaport et al. 2004; Di Salvo et al. in prep.): 

T(t) = M(t) VGMRco --p- (5) 

where M is the mass accretion rate, M the NS mass, // the mag- 
netic dipole moment, Rco the corotation radius (at which the 
disc matter in Keplerian motion rotates with the same angular 
speed of the NS) in cgs units. Since this AMXP is observed to 
spin down, the contribution to the spin up from the accretion 
(first additive term in |5j is neglected. The torque can also be 
written: 

T(t) = 27rlvs(t) (6) 

where I is the NS moment of inertia and vs the spin fre- 
quency derivative. Since vs and Rco are known and, assuming 
I = 10'*^g-cm^, we can obtain the value of the magnetic moment 
and, in turn, an estimate of Bs. Then, through the relation 

B = 3.2 X 10''' ^P Pdip Gauss, (7) 

we can obtain an estimate of the dipolar spin-down Pjip during 
observation (in quiescence), assuming the spin down is governed 
by magnetodipole braking. It leads to Bs = 6.6 • 10^ Gauss and 
Pdip < Ps and then we can conclude that the adopted APtot is 
safely large. Therefore, since the pulse broadening for an error 
equal to 1/20 of the AP interval is less than 0. 1 in phase over our 
observations, we selected 40 trial values of the period, 20 to the 
right and 20 to the left of the nominal value to cover the APtot 
interval calculated above. 

A further check of the validity of the derived Bs (hence of the 
adopted APtot) can be derived from optical observations (Monelli 
et al. 2005; D' Avanzo et al. 2008). In fact, the optical counterpart 
of the XTE J0929-314 companion observed in December 2003 
(very close to our radio observations) with the VLT has a lumi- 
nosity one order of magnitude higher than expected. This lumi- 
nosity excess can be interpreted as the luminosity LpsR isotropi- 
cally irradiated by the rotating magneto-dipole, intercepted, and 
reprocessed by the companion star, as observed e.g. by Burderi 
et al. (2003) and Campana et al. (2004) for SAX J1808.4-3658, 
and by D'Avanzo et al. (2007) for IGR J00291-t-5934. 

The excess luminosity can be written as the fraction f = fc H- 
fn of dipolar radiation intercepted by the companion star and the 
disk (see Burderi et al. 2003): 

924 

Lexc = f ■ LpsR = f ■ = f • 3.85 X lO^^Pl^ erg ■ s"' (8) 

where o) is the rotational frequency of the NS and fc can be 
written as fc - 2n??{\ — cos0)/(47ra^), where a is the the orbital 
separation and 9 the angle subtended by the companion star as 
seen from the central source. If the companion star fills its Roche 
lobe, this can be written as sin Q - RL2/a, where Rl2 is the Roche 

' The expected variation of the Ps during the radio observation is 
neghgible and hence does not affect a correct folding. 
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Fig. 1. Example of a plot resulting from the procedure describee 
in the text. The bottom diagram shows six integrated pulse pro 
files, while the grayscale on the left represents the signal ii 
the 255 subintegrations in which the observation has been sub 
divided. On the right, the parameters used for the folding an 
shown, along with the parameters of the observation. 



Table 2. Parameters for the pulsar candidate with the highest 
S/N shown in Fig.[T] 



Signal to noise ratio, S/N 5.87 

Dispersion Measure, DM (pc cm"^^) 5.5 1 

Orbital period, P„,b (s) 2614.74575 

Projected semi-major axis, a- sin i (It-ms) 6.290(9) 

Spin frequency, vj (Hz) 185.105259 

Ascending node passage. To (MJD) 52405 .48676( 1 ) 




lobe radius of the secondary and RL2/a = 0.49 q^^V[0.6 q^^"^ + 
In (1 + q'^"')] (Eggleton 1983). Assuming a mass ratio of q = 
m2/m2 = 0.02/1.4 and sin i = 1, we obtain fc = 0.003; while fo 
= 0.012 is evaluated adopting a standard Shakura-Sunyaev disk 
model (Shakura & Sunyaev 1973) with rii_io = 2.9, and is given 
by the projected area of the disk as seen by the central source, 
27rRx2H(R) (where R is the disk outer radius and H(R) the disk 
semi-thickness at R), divided by the total area, 47rR^. 

Adopting the value of Lgxc - 3.33 x 10^^ erg ■ s"' reported 
by MonelH et al. (2005) and f ^ 1.15 x 10"^, we can thus derive 
Bs - 7 ■ 10** Gauss, in good agreement with the value estimated 
above. 



3. Results 

Spanning the searched ranges of Porb, Ps and DM we have ob- 
tained 50922 plots, reporting the result from the folding of the 
deorbited and barycentred time series. An example of these plots 
(obtained from the observation at 6.5 GHz) is shown in Fig. 
[1] The bottom diagram shows the integrated pulse profile (6 
phases), while the grayscale on the left represents the signal in 
the 255 subintegrations in which the observation has been subdi- 
vided. A good candidate should display a roughly linear trend in 
the grayscale and a high signal-to-noise ratio pulse profile. On 
the right, the parameters used for the folding of this candidate 
are shown, along with the parameters of the observation. 

A useful diagnostic tool that can be adopted to further eval- 
uate the plausibility of a suspect is shown in Fig. |2] where a 
grayscale of the strength of the signal (with darker points at 
higher S/N), plotted as a function of dispersion measure DM and 
spin frequency vs, should define a clear peak around the correct 
parameters of the putative pulsar The same kind of plot can be 
created at a constant DM with vs and Porb on the axes or at a 
constant vs, tracing the S/N trend at varying DMs and PorbS. 

The plots in Figs. [1] and |2] refer to the suspect with the high- 
est S/N found in our search. Its folding parameters are listed in 
Table |2] The peak in the profile has a 3.4 cr significance that. 



Fig. 2. Signal-to-noise ratio in function of the spin frequency (20 
step above and 20 below the nominal value, corresponding to 0) 
and the dispersion measure (72 steps corresponding to values 
interval from 5.51 to 396.74 cm"^^- pc) at Porb = 2614.74575 s. 
The highest point corresponds to the best S/N, x 6, reported in 
Fig.m 



over the 26568 trial foldings performed on the 6.5 GHz dataset 
(72 DMs X 9 PorbS X 41 Pss), has a probability of not being ran- 
domly generated by noise of ~ lO"*". We also note that, in Fig.|2] 
the decreasing S/N trend is not particularly defined, although a 
maximum is present. Also, that the DM that maximises the S/N 
is close to zero (w.rt an expected value of > 100 pc/cm^*) weak- 
ens the credibility of the signal. Finally, this signal suspect is not 
confirmed in the observations at higher frequency (8.5 GHz), 
elaborated with the same parameters. We can conclude that no 
radio pulsation with the expected periodicity has been found in 
the source XTE J0929-314 during its quiescent phase. 



3. 1 . Flux density upper limits 

A rough estimate of the flux density for a pulsar of period P is 
(e.g. Manchester et al. 1996) 



en„ 



Tsys + T 



sky 



GVNpAtAvMHz VP-W, 



We 



mJy, 



(9) 



where n^- = 6.0 is the threshold (in unit of sigma) for 
having a statistically significant signal, given the number 
of performed foldings; Tjys = 50 K the system noise 
temperature for the observations at ~ 6 GHz and Tjys 
= 25 K for those at ~ 8 GHz (see Parkes website: 



http://www.parkes.atnf.csiro.au/observing/documentationl; Tsky 
is the sky temperature in Kelvin, calculated from that at v = 408 
MHz and considering a scaling with the frequency as v"^^; G is 
the gain of the radio telescope (in K ■ Jy '), with 0.46 for the ob- 
servations at ~ 6 GHz and 0.59 for those at ~ 8 GHz, (see Parkes 
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80 100 

Duty Cycle i%Fs) 



Fig. 3. Flux density upper limits in juJy for XTE J0929-314 ob- 
tained with Eq.|9]for increasing values of the duty cycle. 



website); At is the integration time in seconds, Np the number of 
polarizations (here 2), and Avmhz = 576 MHz is the bandwidth 
in MHz. The term e ~ 1.4 is a factor accounting for the sensitiv- 
ity reduction due to digitisation and other losses. The term We is 
the effective width of the pulse: 



We = -y/w2 + (J36t)^ + 



6tl 



■6tl 



(10) 



Its value depends on the sampling time (6t - 200 fis in this 
case), on the technical characteristics of the receiver (J3 = 2), 
on the broadening of the pulse introduced by both the dispersion 
of the signal in each channel (dtuM ~ 4 -H 9 x lO^^'x DM sec, 
depending on the frequency), on the scattering induced by inho- 
mogeneities in the ISM (^tjcatt ~ 10"^ sec), and on the intrinsic 
pulse width W. The dependence of the minimum flux density 
achievable from the duty cycle, W/P, is shown in Fig. [3] A duty 
cycle of 15% of the spin period is assumed and, with this value, 
the flux density upper limits for XTE J0929-3 14 are: = 68 



at 6410.5 MHz and S 



MAX 



26 juJy at 8453.5 MHz. 



4. Discussion 



In the hypothesis that a radio pulsar was on during our observa- 
tions, we examine several possible explanations for the lack of 
the detection of radio pulsation from our source. 

4.1. The luminosity 

In Fig. |4]the pseudoluminosity L = S d^ in mJy kpc^ (where S 
is the measured flux density and d the distance of the source) 
distribution at 1.4 GHz for a sample of 42 galactic field MSPs 
is showr0. For XTE J0929-314 the upper limits on the pseu- 



^ These values are derived from the ATNF catalogue - 
|http://www.atnf.csiro.au/research/pulsar/psrcat/[ Manchester et al. 
2005 




J 

Log L (8 GHz) Log L (6 GHz) 



Fig. 4. Pseudoluminosity distribution of a sample of 42 field 
MSPs. The arrows indicate the upper limits of the pseudolumi- 
nosity scaled at 1.4 GHz for XTE J0929-314, derived on the 
basis of the minimum flux detectable and assuming a distance of 
6 kpc (Galloway et al. 2002). 



doluminosity has been calculated assuming a distance of 6 kpc 
(Galloway et al. 2002) and our flux density upper limits at 6.4 
and 8.5 GHz scaled to 1.4 GHz assuming a spectral index for 
MSPs of 1.7 (Kramer et al 1998). As a result, about 90% of the 
known MSPs are below the pseudoluminosity upper limits of 
XTEJ0929-314. 

This suggests that XTE J0929-314 in quiescence might not 
be a very bright MSP and that a deeper search should be car- 
ried out to sample lower values of the flux density. However, 
given the distance of XTE J0929-314 (significantly higher than 
the typical distance of the sample of known MSPs), only next- 
generation instruments, like the Square Kilometer Array, will be 
able to discover pulsed radio emission from this source at the 
high frequencies adopted in this work, if its luminosity is at the 
faintest end of the known distribution. 



4.2. The beaming factor 

The emission from a pulsar is strongly anisotropic. This means 
that it does not irradiate in the same way in all directions, but 
only on one usually narrow portion of sky, which can be quanti- 
fied through the so-called beaming factor Therefore, the lack of 
a detection of a radio signal could be due to unfavourable geom- 
etry of the radio emission with respect to the observer. 
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The average value, f(a), of the fraction of the sky swept from 
two conal radio beams of half width a is (Emmering & Chevalier 
1989): 

f(a) - j f (q', 77) sin 77 d/y = (1 - cos Q') + - aj sin Q- (11) 

where 77 is the angle (supposed randomly distributed) between 
the magnetic axis (aligned with the radio beams) and the rota- 
tional axis, and f(a, rf) - cos[max(0, rj-a)] -cos[min(0, 77 + q')]. 

Considering a 10% H- 30% interval of width of the pulse, it 
follows that 9° < Q- < 27° and, therefore, using Eq. [TT] 0.23 < 
f{a) < 0.61. In particular, assuming 15% typical value of width 
of the impulse, the probability that the radio emission cone does 
not intersect our line of sight is 66%, with f(a) » 0.34. 

5. Conclusions 

The aim of this work was to search for radio pulsations from 
the AMXP XTE J0929-314. The detection of radio signals in 
the quiescent phase of this kind of transient systems would be 
ultimate proof of the recycling model, unambiguously establish- 
ing that the AMXPs are the progenitors of radio MSP. No ra- 
dio pulsation has been detected in the analysed data down to a 
limit of 68 yuJy at 6.4 GHz and 26 jiSy at 8.5 GHz. Assuming 
that a radio pulsar was on during the observation, the free-free 
absorption cannot be responsible for the lack of a detection, 
given the relatively high radio frequencies. The beaming factor 
is a viable explanation for that, since the probability of any un- 
favourable geometry of radio emission with respect to the ob- 
server is ~ 60 70%. However, the most likely reason for a 
negative result is that the source has a luminosity lower than our 
limits. In fact, more than 90% of the known MSPs show a lu- 
minosity lower than the upper limits that we have derived for 
XTEJ0929-314. 

To get round to these problems, one can make observations 
at an intermediate frequency between 6.5 and 1.4 GHz. For a 
maximum mass transfer rate in quiescence equal to the outburst 
value, we can estimate that the lower limit in frequency for ob- 
taining Tff ~ 1 is ~ 3 GHz. Calculating Smax for an observation 
with all the same parameters (adopted in the present work) but 
the frequency (now set at 3 GHz), we could sample more than a 
half of the known MSPs luminosity distribution. 
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